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ABSTRACT 


Sitnulation model of various modulation and demodulation schemes used in digital 
communication systems has been implemented using a graphical programming software 
LabVIEW Data generators Ime coders modulators AWGN channel detectors and 
decoders have been modeled for all the modulation schemes simulated A signal fi'om a 
function generator implemented in LabVIEW is the input source signal which has been 
processed and reconstructed 

The effect of channel characteristics m terms of mtersymbol mterference and probability 
of erroi has been highlighted through simulation of zero forcmg preset and adaptive 
equalizers It has been highlighted as to how the Inter symbol interference generated due 
to channel characteristics reduce noise margms and lead to higher errors m detection 

The probability of error plots obtained for various modulation schemes the spectra, 
parameters determined over the link and the results of comparison of various modulated 
carrier signals in terms of power spectra phase transitions or transmission through a 
bandlimited channel match well with the theoretical results 
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CHAPTER 1 


INTRODUCTION 

1 1 LABVIEW & TOOLS OF ANAYLYSIS 

Simulation may be defined as the disciplme whose objeetive is to imitate one or more 
aspects of reality m a way that is as close to that reahty as possible In a simulation 
based approach a computer is used to simulate the waveform or signals that flow 
through the systems Simulation of a commumcation system mvolves generatmg 
sampled values of signals and noise processmg these sampled values through discrete 
time models of functional blocks m commumcation systems and estimatmg 
performance measures 

A modular approach towards a commumcation system would show that it 
comprises of blocks which take mput from precedmg blocks process the signal and 
give their output for further processing These blocks could be a source coder 
modulator, channel filters detectors or decoders If these different blocks could be 
modeled separately and somehow combmed to represent a complete commumcation 
system then it could be classified as a complete simulation The performance measure 
of the communication link simulation is the bit error rate obtained for various 
modulation / demodulation schemes m a AWGN simulation channel For baseband 
schemes the Eye diagram is taken as a subjective measure of performance 

Digital simulation can provide a useful and effective adjunct to analytical 
performance evaluation or direct hardware evaluation of commumcation systems It is 
a computationally demanding task due to the large number of repetitive signal 
processmg operations that must be performed m order to obtam a statistically valid 
measure of system performance Various computer languages have been used in the 
past like Simula, FORTRAN Pascal to simulate commumcation systems, usmgWICS 
or SKY VORTEX processors Due to the technological advancements, today a 
personal computer has the capabihty of such analysis and simulation usmg suitable 
software packages The software package used here for simulation is LabVIEW 
(Laboratory Virtual Instrument Engmeermg Workbench) LabVIEW is a program 
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development application much like C or BASIC havmg one basic difference that other 
programming systems use text based languages to create lines of code while 
LabVIEW uses a graphical progiammmg language G to create programs m block 
diagram form 

The mam aim of the thesis has been to build mdependent menu driven blocks 
where the user can decide the nature of the commumcation link he wants to simulate 
by selectmg appropriate modules At each stage the processed output is available for 
viewmg graphically so as to understand better how the waveforms behave at various 
stages of a link The complete Passband and Baseband simulations can be run from 
two Panels shown m Figure 1 1 Each control knob when imtiated gives further 
options for the user to select from 
The AIM of the thesis can be defined as 

1 To convert block concept level commumcation Imk mto a modular 
menu driven utihty which is user mteractive 

2 To construct modules of various modulation demodulation schemes 

3 To simulate channel equalization and its effects 

4 To assess performance of these detection schemes usmg the BER 
plots 

5 To give the user a stage wise graphical representation of processmg 
over the link 

To understand the advantages of LabVIEW let us first understand how the 
programmmg is done There are three basic elements necessary to build any source 
code (i) Tenninals (ii) Nodes and (m) Wires The mput parameters are operated upon 
by the flinction chosen and the output of this function block is the desired result The 
source code which is the Block diagram is hidden from the user and he has access only 
to the front panel where he can specify the mputs and see the results The followmg 
can be classified as the mam advantages of this software package 

(1) Memory Usage One of the mam challenges of conventional language is the 
allocation and de allocation of memory LabVIEW a data flow paradigm of C 
removes this difficulty as you do not need to allocate variables nor allot values to and 
from them Instead you can create a diagram with connections representmg the 
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transition of data Functions that generate data take care of allocating the stoiage for 
the data When data is no longer used the associated memory is deallocated 

(2) Debugging Techniques & Performance Profiling LabVIEW has strong 
debugging tools like probes (to see values of data at intermediate points) breakpoints 
step runmng of program, highlight option which runs a VI in a slow mode highlightmg 
all data values An Error VI mdicates the various programmmg errors due to which 
the program is nonfunctional It also has performance profiling which is an mteractive 
tabular display to view the performance data of sub Vis when called firom a specific 
VI 

(3) Graphical Displays and Data Acquisition LabVIEW (version 5 0 onwards) 
provides six types of graphs and two types of chart displays These graphs include the 
3D surface graphs parametric graphs XY plots waveform graphs and waveform 
charts Data at any intermediate pomt m the hnk could be easily given for plottmg 
Seeing output parameters visually without any complex source codmg is a big asset 
associated with LabVIEW System parameters can easily be analyzed at intermediate 
steps if they are graplncally represented 

Data Acquisition (DAQ) card facilitates a high level mterface to DAQ devices and 
signal conditionmg hardware NI DAQ driver software interfaces the LabVIEW Vis 
with DAQ hardware devices and helps m developmg complete mstrumentation 
acquisition and control application Analogue or digital signals can be given as mput 
parameters to LabVIEW for processmg and also the output parameters through the 
mterface to hardware devices Real tune signals can be processed and analyzed usmg 
these DAQ mterfaces 

Another advantage of LabVIEW over other text based languages is that 
multiple processes can be mitiated at one tune This feature not only helps m 
synclironization but also helps m reducmg the total processor time to complete a 
computation 

(4) Frequently Used Functions 

(a) String to decimal String function used m Ime code implementations for 
searchmg certam patterns of bit sequences 
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(b) Index Array If the output vector is one dimensional mdex array gives the 
parameter value at specified index 

(c) Subset Array Returns a portion of the mput vector 

(d) Sine wave Generates sine waves of specified amphtude phase and frequency 

(e) Integral X(t) Performs the discrete integration 

F(t)= Jf(t)dt 

Integral x = y, = )dt 

o j=0 

Convolution h(t)= x(t)*y(t) where * stands for convolution 

00 

= Jx(T:)y(t-T)dx 
-00 

n-1 

h. = Z^kYi-k 1 = 0 1 2 n 1 

k=0 


00 

(g) Cross correlation R^y (t) = x(t)* y(t) = J x(T)y(t + t)dT 

—00 

where * stands for convolution 
n-1 

hj= Ixkyj+k 

k=0 


Rxy, == h, („ ,) for i=0 1 2 (n -1) 

(h) FFT Computes real FFT or real DIT of mput sequence 

(i) Amplitude and Phase spectrum Computes smgle sided scaled amplitude 
spectrum magmtude and phase of a real time domain signal 


()) Butterworth filter Can be used as a low pass/band pass filter Samplmg 
frequency low and high cut off frequencies and the order have to be specified 


(k) erfc(x) 


00 


2 0 
j exp(-F^)dt 

Vxx 


4 



(l) Standard deviation = —'^(x 

I n-1 

where y= — n = number of elements 

^^ 1=0 

(m) Spline interpolation used for interpolating discrete pomts for graphical 
representation 

(n) Solve linear equations solves a real linear system AX=Y 

where A X and Y are vectors 

(o) Create special matrix To generate nx by ny toephtz matrix 

(p) Analogue Input acquire waveforms used in DAQ programs Acquires data 
from the specified channels and samples the channels at the scan rate 

12 GENERAL DESCRIPTION OF COMMUNICATION 
SYSTEMS 

The schematic of a basic commumcation system shown m Figure 1 0 have been 
coded using LabVIEW Detailed and software implementation are illustrated m the 
subsequent chapters A brief review of some basic concepts are appended below 



Figure 1 0 Basic communication System Block 











Sampling 


The sampling theorem states that a band hnuted signal having no spectral 
components above fm hertz can be determined uniquely by values sampled at uniform 
mtervals of Ts seconds where 


Ts<l/(2fn) 

The upper limit on Ts can be expressed m terms of the sampling rate denoted 

fs-l/Ts 

T^fm 

The sampling rate fs = 2fm is also called the Nyquist rate [10] 

Analogue to Digital conversion 

The process of digitizmg the original analogue signal is called Quantisation It 
consists of mapping the amphtudes of the signal mto a number of discrete amplitude 
levels This process results m loss of information smce it is not possible to recover the 
origmal analogue signal [1 1] 

Orthogonality Principle 

Signals that satisfy the following condition are orthogonal to one another 

T 

\s,{t)Sjit)dt = ES,j 

0 

where 8,j = l 

= 0 i=j 

An N dimensional orthogonal space can be defined as a set of N Imearly independent 
functions (jijCt) called the basis function Each (|)j(t) function of the set of basis function 
to satisfy the principal requirement for orthogonahty must be independent of the other 
members of the set and it must not mterfere with the detection process of the other 
members of the set [13] 

Probability of Error in AWGN Channel 

s(t)= Acosoot 
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IZE 

= - — cosaf 
V T 

where P repiesents the average Power in watts and E the energy per bit [22] 

Eb _ ST _ S _ SW Sj"W^ 

No No RNo RNoW“nIr, 

Eb = signal energy per bit 
No = noise spectral density 

S = signal power 
T = bit duration 

R = Rate 

W = bandwidth 

A = Amplitude 

CO = 2ix(lTequency) 

Correlation Receivers 

A matched filter and a product mtegiator both perform the same task v 
T 

Z(T)= J r(T)s(T)dT where x is the time difference 
0 

The product mtegration of the received signal r(t) with a replica of the transmitted 
signal s(t) IS known as correlation of r(t) with s(t) At time t = T the output of the 
correlator and matched filter is identical 
Response of Correlators to White noise 


where No/2 is the power spectral density and a/ is the noise variance [22] 

ISl and Equalization 

A band limited channel spreads a pulse waveform and causes pulse overlap 
This overlappmg of pulses is called Intersymbol interference (ISI) An equahzer filter is 
used to compensate for the distortion caused due to the channel The minimum 
bandwidth needed to detect Rs symbols/sec without ISI is IV2 Hz If the rate 
mcreases beyond 2B where B is the bandwidth it effects the discrete samples given to 
the sampler for detection ISI reduces the noise margins and can lead to higher errors 
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The process of correcting the channel distortion is called equalization A (2N+1) 
tap transversal filter with optimal tap coefficients determined by the peak distortion or 
mean square distortion criterion forces the equalizer output to zero at N pomts on 
either side of the desired peak The weights m MSE multiplied with the equalizer 
mput adds up to the polar response of the symbol ( 1 or — 1) bemg sent In both cases 
N snnultaneous equations are solved [21] 

13 PROBLEM FORMULATION 

The thesis mvolved simulation of the generic communication system model 
usmg LabVIEW package and the evaluation of system performance with various 
modulation demodulation schemes The thesis covers the following 


(l) Passband Systems 

(a) Generation of an analogue source 

(i) Generation of smusoid signal at various frequencies amplitude and phase 

(ii) Generation of random signal 

(ill) Signal taken from function generator using DAQ utility 

(b) Conversion of analogue signal to digital using linear and non Imear 
quantisation techniques 

(i) 4 level Imear quantiser 

(u) 7 level Imear quantiser 

(m) A law non Imear quantiser 

(iv) [I law non Imear quantiser 

(c ) Lines codes 

(i) Alternate Mark Inversion coder (AMI) 

(ii) m Bmary n Binary coder (mB nB) 

(ill) Bipolar Six Zero Substitution coder (B6ZS) 

(iv) Biphase coder 

(v) 4 Bmary 3 Ternary coder (4B3T) 

(d) Simulation of various digital modulation schemes 

(i) Amphtude shift keymg (ASK) 

(u) Coherent and non coherent frequency shift keymg (CFSK Sc NCFSK) 

(m) Bmary phase shift keymg (BPSK) 

(iv) Quardrature phase shift keymg (QPSK) 

(v) Mmimum shift keymg (MSK) 

(vi) 16 QAM and 64 QAM 

(vii) 3PSKand3FSK 

(vm) Offset quardrature phase shift keymg (OQPSK) 


8 



(e) Simulation and generation of AWGN samples 
(:0 Simulation of demodulators (of all modulators) 

(g) Reconstruction of input signal after decoders (using different coders and 
quantisers) 

(h) Evaluatmg system performance for different SNR m terms of probability of 
error 

(i) Determimng spectrums and parameters over the link 

(I) Evaluation of the effect of bandpass filters on QPSK and OQPSK signals 

(II) Baseband Systems 

(a) Generation of bmary source sequence 

(b) Perionnance of zero forcing equalizer using different channel impulse 
response 

(c ) Performance of a preset equalizer and decision directed Imear equalizers 

(d) BER evaluation of decision directed Imear equalizer 

(e) Channel tracking and errors m decision directed Imear equahzer 

(f) Generation of Eye diagrams before and after equalization 

The above operations mclude all functions of transmitters AWGN channel and 
receivers A known bit sequence is transmitted duly corrupted and received The 
comparison m terms of probabihty of error of the transmitted and the received bits 
gives the system performance 

1 4 Organization of Thesis 

After the mtroductory chapter deahng with advantages of LabVIEW, basic 
description of commumcation system model and underlymg prmciples Chapter 2 
explams the various Passband modulation schemes data generation line codmg 
quantisation, decoders and demodulators Chapter 3 deals with Baseband schemes and 
Equalization Chapter 4 brmgs out the various results like the probabihty of error m 
various modulation schemes errors caused due to ISI errors caused due to bandpass 
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filtering eye diagrams and certain parameters determined over the link Chapter 5 
deals with the concluding remarks and suggestions for fiiture work 

Some of the unportant programs used in the study are appended m the 
Appendix The full package consisting of 300 files (65 MB) is submitted on a CD 
along with the thesis report 
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CHAPTER 2 


PASSBAND MODULATION AND DEMODULATION 

SCHEMES 

2 10 DATA GENERATION AND CODING 

The generation of analogue waveform is simulated m three different ways 
(i ) Smusoidal waveform generation usmg a function generator usmg function vi 
(u ) Takmg an analogue mput signal from a function generator (external source) using 
Analogue Input vi 

(ill ) Random signal generation usmg random vi 

Figure 2 1 shows the front panel of the function generator The parameters available to 
the users are the frequency, phase and amphtude The generated waveform (sme 
sawtooth or square) can be visually seen on the graphical display The frequency is 
defined by the samplmg rate specified Upto fjl Hz of frequencies can be generated on 
the function generator where f, is the samplmg frequency The frequency m cycles/sec 
divided by the samplmg rate gives the normalized frequency[2] 

cycles / s _ cycles 
samples / s sample 

The reciprocal of the normalized frequency gives the no of samples requued to 
generate one cycle of the sme wave By controllmg the fine and course knobs we can 
generate frequencies upto ij 2 A digital display facihtates the user to see what frequency 
has been generated much like checkmg the frequency on a CRO The amplitude and phase 
can also be adjusted DC offset or noise can also be added to the generated waveform 

The random number generator function is utilized to generate a random waveform 
pattern (uniformly distributed) The user specifies the range that is the maximum and 
minimum values of the voltage 
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The third option is to take m a signal from an external source a function 
generator, llirough the DAQ interface The scan rate specified in the program specifies 
upto to what value the frequencies can be generated The scan rate is a variable 
parameter Acquire waveform vi using the mtermediate level DAQ vis is used to acquire 
waveforms for further processmg The select mputvi is a menu driven utility which 
selects the desired input from the above options 

211 Quantisation 

The conversion of simulated analogue waveform to discrete values is done via the 


quantisation vi 

These are classified mto four quantisation processes 

(1) 

16 level Imear quantisation 

(2) 

128 Level Imear quantisation 

(3) 

A law nonhnear quantisation 

(4) 

u law nonhnear quantisation 


Implementation of 128 level linear quantisation program (Program #1) 

1 Determme the range (maximum voltage minimum voltage) 

2 Divide range by 128 

3 Run iteration x2'’ 7 times givmg output (2^ 2^,2'* 2^ 2^ 2' 2°) 

4 Multiply the above array o/p with the o/p of (range/1 28) We get the X* Y array 

5 Classify mput as + or and give mput to appropriate case statements 

6 Compare first value of X*Y array with the mput, if the mput value is greater, subtract 
the two and now compare to the second value of the X*Y array Assign Boolean 
True if \alue is smaller, or False if greater Use slnft register to store added values 

7 Carry out above operation for all 7 bits and get an array of Boolean true and false 

8 Convert Boolean to 0,1 (1= true 0= false) 

9 Each mput voltage sample is converted to its correspondmg seven bits 

10 For negative mput carry out above conversion to Boolean m the exact manner 
described usmg the case statement for the -ve mputs 
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The above is the 2 s compliment method of Quantisation To illustrate numerically let 
the input voltage value be +4V (The peak mput voltage is +5 V and the number of levels 

IS 16 ) 

10/16=0 625 

[X*Y]=[5 2 5 1 25 0 625] 

Compare 4 with 5 0 assign False is 4 >5 

Compare 4 with 2 5 assign True is 4 >25 

Compare (4 2 5) with 1 25 assign True is 1 5 >1 25 

Compare (4 2 5 1 25) with 0 625 assign False is 0 25 >0 625 

So we get for an input voltage of 4V a bmary of 01 10 (Refer to Table 2 1) For some 
range of mput voltage values a bmaiy number is assigned This range of mput voltage is 

P 2V 

called the step size [q = — = — ] Here V=5 volts and M is 16 The voltage values are 

M M 

assigned m steps of 0 625V[10] 

NonLinear Quantisation 

Tapermg of the quantisation level leads to better performance as compared to the equal 
spacmg quantisation The mput signal is non hnearly compressed m amphtude to force all 
signals to he withm a specified range This compression is typically of logarithmic form 
The fj. Law and A Law are two such classifications [16] 


p Law 

, , ln(l + plxl/V) , , 
y(x) = — _LJ_sgn(x) 
in(l + p) 


xUl 


A law 

Ax 

I + lnA 



_ 1 _ 

A 
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, _ 1 + lnAx ^ ^ 1 1 , 

y(x) = — _sgn(x) 

1 + lnA A ' 

Typical value of A being 100 and p is 255 Where y(x) is defined as the quantised level 

assigned for x the normalized mput signal 

VI for p law quantisation is described below 

1 Sign magmtude format is used where the first bit gives the polarity The remaining 7 
give the magmtude (Refer to Table 2 2) 

2 Assign 1 if polarity is positive Assign 0 if polarity is negative 

3 Next find S the segment identifier If input voltage between the range 0 to 3 1 it is 
mapped to segment 1 

x<64 2“ 33 X = sample value 

a=0 1 7 

4 Next determme residue r that is the difference between input amplitude and the lower 
end pomt of the segment 

r=x s=0 

r = x(32 2“ 33) s-12 7 

5 Next determme b the quantisation interval contammg r 

r < 2b+l 

r<(2^"')(b+l) s=12 7 

6 On gettmg polarity S and b assemble the 8 bits Sample mput has been mapped to 
bmary 8 bits 

7 Invert bits for transmission The range of encodmg Table is upto 8159 If mput 
voltage value exceed tins limit an error signal is generated Illustratmg the above 
numerically 


Input voltage = 24 volts 
+ve Assign 


1 
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Segment 1 Assign 000 

bis 12 Assign 1100 

Biiidiy 10001100 

A law follows similar algorithm [16] 

212 LINE CODE 

Coding is the transformation of the source bits to the transmitted data symbol The 
coders add redundancy to the transmitted data symbols and at the receiver end checks 
are made to see wheather these code constraints are preserved or not Some of the 
advantages of coding are 

1 To place spectral nulls at dc thereby enabling baseband transmission through 
a c coupled channels 

2 To add redundancy and detect errors 

3 Line codes allow concentration of signal power at frequencies near dc 

4 Lme codes have effective control over the power spectrum 

5 Lme codes assist in extractmg timmg or clock information from the signal 
The followmg codes have been implemented which form the codmg module These menu 
driven options are 

1 Biphase coder 

2 Alternate Mark Inversion coder (AMI) 

3 Bipolar Six Zero Substitution coder (B6ZS) 

4 m Bmary n Binary coder (mBnB) 

5 4 Bmary 3 Ternary coder(4B3T) 

Of these mBnB codmg and its miplementation is discussed below 
mBnB coder 


DISPARITY AND DIGITAL SUM 

The disparity of a digit ( 0 1,2 r 1) is defined as[15] 

d, = a, (I) 
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(II) 


where ai is the ith digit, is the complement of it, and r is the radix 
^ = (r I)-a, 

The disparity of an n digit code word is given by 

n n 

Edi =2 -n(r-l) 

1=1 1=1 

n 

= 2 2] £ 1 , n (for binary transmission) 

;=1 

= (No of 1 s) - (No of 0 s) 

The digital sum ds, of an n digit code word is given by 

/=i 

_ J -0 5 for a zero 
I 0 5 for a one 

In the mBnB code the bmary sequence which is to be transmitted is first broken mto m 
bits and then converted to n bits and transmitted (where n > m) The number of 
combmations of each code havmg a particular disparity (assummg n = m +1 and m is odd) 
is given by 

No of words havmg zero disparity = C(m+iy 2 

No of words havmg Non zero disparity = C (m+ lya 

A source word is translated to zero disparity or a pan: of words havmg alternate disparity 

3B4B code 

3 bit mput data is converted to 4 bit output data 
RDS max = 1 

RDS mm = 1 

where RDS (Runnmg digital sum) at the k"’ bit is defined as 
k 

RDS(k) = la, 

1=— oo 

DS V (digital sum variation) = RDSmax — RDSmm 

-2 
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the coding done should give the minimum DSV 


IMPLLMEN lAllON (Program U 2) 

1 Initiate all mput variables to imtial value 

2 Reshape mput array from 2 dimensional to one dimensional array of 3 bit serial input 
using the reshape array function 

3 Refer to Table 2 3 Map 3 bit serial mput dnectly to 4 bit bmary codeword if the code 
and its complimentary are same 

4 If not same, calculate disparity If DSV > 2 then select ve disparity code and if 
disparity < 2 then select +ve disparity code This is implemented by using While Loop 
for disparity calculation and then usmg appropriate case statements that do the 3B to 
4B mappmg 

5 Repeat step 2 tlirough 4 for all input bits and reconstruct the mapped one dmiensional 
array of the converted 4 bits 

6 Give this output to chart for plotting 

BIPOLAR SIX ZERO SUBSTITUTION (B6ZS) 

To modify the Imear AMI Ime code to a nonlmear code, B6ZS block Ime code is 
used This class of code modifies the AMI code by performmg a substitution for a block 
of SIX consecutive zeros For each block of six zeros the code word BOVBOV is 
transmitted where B stands for bmary and V for violation If the RDS at the start of the 
block is zero then, +0+ 0 0 would be substituted for the six zeros, otherwise -0 +0+ 

4 BINARY 3 TERNARY 

In this block code bmary words are mapped mto ternary words, with the use of 
feedback to choose the appropriate mode at each block so as to minimize the RDS 4B3T 
IS a bimode code where mode A is chosen whenever the RDS at the begmnmg of the 4 bit 
block IS m range 3<RDS:S 1 and mode B is chosen when the RDS is m range 0<RDS:^ 
Table 2 4 has been used to implement 4B3T coder 
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2 2 MODULATION SCHEMES 


The followmg modulation schemes have been mcorporated usmg Lab VIEW A 
buef theoretical preview will follow the actual implementation 

1 Amplitude shift keying (ASK) 

2 Coherent and Noncoherent frequency shift keymg (FSK& NCFSK) 

3 Bmary phase shift keymg (BPSK) 

4 Quararature phase shift keymg (QPSK) 

5 Mmimuin shift keymg (MSK) 

6 1 6 Quardrature Amplitude Modulation (1 6 QAM) 

7 64 QAM 

8 3 Phase Shift Keymg (PSK) 

9 3 Frequency Shift Keymg(FSK) 

1 0 Offset Quardrature Phase Shift Keymg (OQPSK) 

The parameters used m the theoretical expressions for the followmg sections are defined 
as 

Eb = Signal energy per bit 

Tb = Bit duration 

CO = 27tf where f is frequency 

^(t) = Basis function 

h = deviation ratio 

A = Amplitude 

0(t) = phase 
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2 2 1 Binary Phase shift keying 



A 


Fig 2 2 BPSK TRANSMITTER 


The general analytical expression for PSK is [14] 

cos[cDot + 0i(t)] Ost:<Tb 


Pb 


s. (t) = "V Tb 


1 = 01 


0i(t) = 


2ni 


M 

= 7tl 


M = 2 



where (t)i(t) = J— cos(27rfct) 

Tb 

The data signal shills the phase of the waveform from one of the two states that 
IS, either zero or % The signals are said to be antipodal (Refer Figure 2 2) 


Implementation 

1 Generate carriers that are offset by 180 degrees using the PSK carrier vi 

2 Depending on the input bit, shift the phase of the carrier from 0 to 180 degrees 

3 Plot output on chart Give output to global modulated 0/P’ 


2 22 Non Coherent FSK 

f A cos(2t^ /) O^t^T, 

1 0 elsewhere 

f equals ft and f 2 
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f; = n,/Tb where n, is an integer and Tb is the bit duration The two frequencies are 
orthogonal (They satisfy the othogonahty principle) fr represents a 0 and f 2 represents a 
binary 1 [11] 

Implementation 

1 Select fc and using h = 1 ( h is the deviation ratio)and R (Bit rate) generate fr and fi 

2 Using nfsk carr vi generate carriers at frequencies fr and fz 

3 Modulate carriers with the incoming bits and store in global "non coh o/p ' for further 
processing Give to chart to plot output 

2 2 3 QUARDRIPHASE SHIFT KEYING 

The mformation carried by the transmitted signal is carried in the phase[13] 


s, (t) = 


Acos(2Tcfct + (2i - 1)^) 

' 0 


0<t<T 

elsewhere 


Si (t)= ^^cos[(2i-l)-^]cos(27cfct)-y— sm[(2i-l)-j]sin(27rfct) 

where i = 1 2 3,4 E = symbol energy T = symbol duration 
Orthonormal basis functions are 

<*^i(0 = 'j^cos(2;r/c0 
f^2(0 = J~sm(2;fcO 


21 




Fig 2 3 QPSK TRANSMITTER 


Signal vectors associated with the four message points are 


Si(t) = 


VEcos((2i-1)-) 

4 

-VEsm((2i-l)-) 

4 


where 1 = 12 3 4 


(I) 


Implementation 

1 Imtiahze all mput variables to default values 

2 Generate (jii (t) and (jia (t) 90“ offset from one another at fc 

3 Generate signal space diagram by varymg i m the above equation (I) and give to an 
XY graph for plottmg 

4 De mterleave the mput bits mto odd and even bits 

5 Generate (jii (t) or <1)2(0 over two bit interval (2T) for the Inphase and Quirdrature 
channel 

6 Add the Inphase and Quardrature channel signalse, the resultant is the QPSK signal 
waveform 

7 Plot the Inphase, Quadrature and the resultant QPSK waveform for the complete set 
of mconimg bits 
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2 2 4 MINIMUM SHIFT KEYING 

In QPSK due to the alignment of the mphase and quardrature (ai (t) and aj (t)) 
streams the earner phase can change over eveiy 2T A change m both components 
simultaneously can result m a phase shift of 180° This phase shift m band limited 
operations results in the envelope gomg to zero The restoration of which leads to 
restoration of all frequency sidelobes back to their origmal level To circumvent this 
problem the modulation scheme minimum shift keymg was developed , where the phase 
of the outgomg wave from the transmitter has no transitions it is contmuous m phase [24] 

s(t)= ^■^^[cos(27cf^t4■6(0)] for symbol 1 

V m 

s(l) “ ^[cos(27tf2t + 0(O)] for symbol 0 0<t<Ti, 

V 'll) 

where 


0(t)=0(0)±^t 


fo + 



fc 



0:St£Tb 


or 

£=i(fi+f2) 

h-Tb(fl-f2) 

where h is the deviation ratio (h = ‘A for MSK) 


The orthonormal basis functions cpi (t) and 92 (t) are defined as the smusoidally modulated 
quardrature carriers 

(p , (t) = cos(^ 0 cos(2jf ct) 0 < f < 7 ; 

■\T b 
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(p2 (t) - smilftfct) 0<t<Tt 

S(t) = S| 9i ( (t) + S2 (t) 92 (t) 

n 

Where si= ji(()<Pi(Odt 

T 

2Tb 

S2= \s{t)(p^{t)dt 


Signal Space Characterization of MSK 


Transmitted Bmary Symbol 

Phase States 

Coordmates 


0(0) 

0(Tb) 

Si 

S2 

1 

0 

+Ttl2 

+VEb 

"V^b 

0 

n 

+n/2 



1 

n 

n/2 



0 

0 

n/2 






Table 2 5 

The above table is used m selecting the phase and the polarity of the scaled time function 
Si 9i(t) and S2 92 (t) Si 9i(t) and S2 92 (t) are offset by Tb The MSK waveform is 
generated by addmg these two waveforms This resultant waveform is contmuous in 
phase 

Implementation (Program # 3) 

1 MSK parameter vi determines the fl and f2 carrier frequencies The input parameters 
required are fc, R(ratc) and h (h=l/2) 

2 MSK trelhs vi generates the phase trelhs usmg the mput bits A nil transition is 
defined for bit 1 and nil transition for a bit 0 The trellis is plotted on a XY graph 

3 Generate MSK orthonormal basis functions 91 (t) and 92 (t) usmg the MSK carrier vi 
generatmg carriers at 180° oJTset for 9i(t) and 92 (t) 
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4 De interleave the input bit stream into odd and even bits (Inphase and Quardrature) 

5 Determine polarity of Si for the mphase bits using the signal space Table 2 5 usmg the 

incoming bit and 0(0) 

6 Determine polarity of S 2 for the quardrature bits using the ineommg bit and 0(rb) 

7 On getting the polarity (positive or negative) of si and S2 generate si(pi(t) and S2 92 (t) 
both over 2Tb, S2 92 (t) bemg offset by fb 

8 Add the Inphase and Quradrature channel waveforms The resultant is a MSK 


waveform 

Thus MSK modulation determines the polarity of the carrier by seemg the previous phase 
The four possible cases can be 


1 

0(0) = 0 

0('Ib) = 7c/2 

BIT 1 

2 

11 

o' 

0 

0(Tb) = 71/2 

BITO 

3 

0(0) = % 

0(Tb)= 71/2 

BITl 

4 

0(0) =0 

0(Tb)= 71/2 

BITO 


▼ 


a. (t) 



MSK 

Signal 
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225 DIFFERENTIAL PSK 


DPSK consists of two basic operations differential encoding of the binary input 
sequence and phase sluft keying I o send a symbol 0 phase advance the current signal 
wweform by 180 degrees and to send a symbol 1 we leave the phase unchanged The 
relative phase difference between two successive bit intervals is stored inDPSK[17] 

To encode the incoming bit stream a reference symbol is the arbitrary first bit Let 
bK be the mconimg bit sequence then the differentially encoded bits are the complement of 
the module 2 summation of the bk and dk i (Refer to Figure2 5) The differentially 
encoded bits modulate a carrier wave fc to produce the DPSK signal 



Fig 2 5 DPSK TRANSMITTER 


DPSK 

Signal 


Implementation 

1 The mitial reference bit is taken as 1 


2 The input bit stream bk is converted to Boolean and XOR operation is carried out with 


the precedmg bit i e , dk i to generate dk (refer Table 2 6) 
3 DPSK carrier is generated usmg DPSK carr vi 


bk 


1 

0 

0 

1 

0 

0 

1 

dk 1 


1 

1 

0 

1 

1 

0 

1 

dk 

1 

1 

0 

1 

1 

0 

1 

1 

I ransmiKed 

phase 

0 

0 

7t 

0 

0 

7U 

0 

0 


Table 2 6 
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4 The encoded bits modulate the carrier i e a change from 1 to 0 the phase changes by 
% and a change from 0 to 1 the phase sent out is 0 For binary 1 the phase is 0 and for 
binary 0 it shifts to n 

5 The DPSK waveform is plotted on chart and output is given to global variable 

modulated 0/P ' 

226 64 QAM TRANSMITTER 

In the M ary Quardrature amphtude modulation scheme the carrier experiences 

amplitude as well as phase modulation The general form of M ary QAM is defined by 
Si(t) = Aaj cos(27rf(;.t) + Ab^ sin(27i:fct) 0< t :^T 

a, = ±1 ±3 Vm 1 b, = ±l ±3 Vm 1 

where a, and b, are chosen in accordance with the location of the message point Each 

constellation point m the square constellation is represented by a unique 6-bit symbol 

which IS Gray coded to mmirmze the error probabihty[19] 

Implementation (program #4) 

1 Generate the square constellation signal space diagram for 64 QAM usmg space vi 

2 6 bits which form a symbol which are given as mput to the For Loop The For Loop 
does the followmg repetitive processmg on all the 6 bit mput arrays 

3 The 6 bits are further spht into two groups of 3 bits by the Array subset function 

4 The decimal to strmg conversion is used to convert to strmgs 

5 Refer to square 64 QAM phasor constellation diagram (Figure 2 6) Group I of 3 bits 
IS mapped to the first 3 bits and Group II to the last 3 bits of the 6 bits assigned to 
the message points of the constellation These bits are gray coded 

6 On domg the mappmg and knowmg where the message point is determme a, and b , 
which would be 1, 3 5 7 1 3 5 or 7 on the real and imagmary axis 

7 Generate QAM carriers Acos(2Tcfct) and Asin(27i:ft) usmg QAM carrier vi 

8 Multiply a. with Acos(2nfct) and b, with Asm(2nfct) and add the two The resultant is 
the QAM waveform 


27 




9 Plot QAM waveform on chart and give modulation 0/P to global variable 
modulated 0/P 

NOTE 16 Q/^M is also implemented m a similar manner 

2 3 DEMODULATION TECHNIQUES 

The next block considered after the transmitter is the receiver The channel is 
assumed to be an ideal AWGN Channel which adds gaussian noise to the transmitted 
signal We will break up the receiver into two parts First we wdl receive input back to 
back adding no noise Next we will see the performance of these receivers under varying 
AWGN variance Tins latter part will be covered under the result section The receivers 
of all modulation schemes discussed in section 2 2 are developed using either the 
correletoi detector or matched filter Amongst these illustration and implementation of 
the followmg receivers is given 

1 Bmary phase shift keymg 

2 Quardratuie phase sluft keymg 

3 Minimum sluft keying 

4 64 Quardrature amplitude modulation 

5 Non coherent firequency shift keymg 

(Proper camerj phase and bit tuning synchromzation are assumed m the receivers for the 
various coherent schemes ) 

2 3 1 Binary Phase Shift Keying 

A pan of signal si (t) and S 2 (t) is used to represent a bmary 1 and bmary 0 

51 (t) = cos(2nfct) 

52 (t) = ^ cos(27tfct + 7t) = - cos(27rfct) 

V Tb V Tb 

The basis function 


(p, (t)= jAcos(27rfct) 



The received signal is 

r(t) = s(t) + n(t) where s(t) is the signal and n(t) is AWGN 

The input signal and the locally generated cpi (t) are available at the product 
modulator The output of the corrcletor (Integrate and dump) is used for comparison with 
a tlireshold level If the correlator output is greater then threshold decide for Binary 1 
otherwise a Bmary 0 [22] 



Implementation 

1 The ’modulated o/p global variable is the mput signal to the receiver The For Loop 
IS used to sequentially give the mput pertaimng to one bit at a time for further 
processmg 

2 The locally generated reference signal is polymorphically multiplied to this mput 
waveform 

3 The resultant is given to the mtegration function 

4 The value after the correlator is sampled at T and is compared to a threshold value 0 
If greater than 0 the case structure gives a 1 If less than 0 than the case structure 
gives a 0 as output 

5 The O/P bits are the same as the quantised bmary mput to the transmitter 

6 Ihe output bits arc given to the "rcc O/P' global variable lor further processing to 
follow m the decoder 
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232 QPSK RECEIVER 

The QPSK signal as defined in section 2 2 is given by the equivalent form in the 
interval 0^t<T by 

r 

K 


2E 

s.(t)= J-^co$ 


( 2 i-l)- 


[2E 

. _ ^ . 7t 

1 — sm 

( 2 i-l)- 

^ T 

L 4 J 


sin( 2 Kfct) 


and the two basis functions 


91 (0 = 

92 (t) = 



cos(27tfct) 

smilnfct) 


0<t<T 

o</<r 


The incoming bit stream at the transmitter was partitioned mto the Inphase and 
Quardrature channels The QPSK receiver consists of two correlator channels with the 
locally generated 91 (t) and 92 (t) basis functions The output of the correlator of each 
channel is compared to thresholds and the o/p bits are decided based on this comparison 
These output bits are then multiplexed to give the origmal input bmary sequcnce[13] 


r(t) 



Fig 2 8 QPSK REC 


Implementation 

1 Convert QPSK signal to a 2 dimensional array form and give to For Loop to process 
one T sec waveform at one time 


Binary 

Sequence 

-> 
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2 Generate (pi (t) and (p 2 (t) for the Inphase and Quardrature channels 

3 Polymorphically multiply the received signal with the cpi (t) and (pa (t) 

4 Integrate the resultant one dimensional array in both channels 

5 Index the output of the mtegrator at T and discard all other values 

6 Give output of the For loop to mdicators labeled mphase and quardrature 

7 Give the mdicator 0/P through a local variable to a For Loop havmg a case structure 
which gives a Bmary 1 if mput is greater than threshold otherwise 0 The threshold is 
set to zero 

8 The binary output is stored m mdicators labeled Even and Odd 

9 Use the mterleave function usmg Local variables Even and Odd multiplex to get the 
desiied output Tins is the detected bmary sequence Give this output to chart for 
plottmg 


233 64 QAM RECIEVER 

The transmitted waveform is given by [19] 

Si(t) = Aai cos(27i:fct) + Abj sm(2Tcfct) 0<t :^T 

The amphtude 7d, 5d, 3d, d, d 3d 5d and -7d of the I and Q signals are 
assigned 3 bit gray codes 011 010 000,101 100 1 10 and 111 respectively The three I 
and Q bits are denoted by ii I2 13 qi q2 qa respectively These mterleaved bits are 
demodulated usmg the decision boundaries The recovery of the bits is based on 
I Q>0 ii qi =0 

I Q < 0 ii qi = 1 


I Q>4d 

12 

q2 = 1 

4d^I Q<4d 

I 2 , q2 = 0 

4d>I Q 

12 

q2 = 1 

I Q>6d 

I 3 

q3 = 1 

2d<I Q<6d 

I 3 

qs =0 

2d<I, Q<2d 

13 

q3 = 1 

6d<I, Q< 2d 

13 

q3 =0 

6d>I Q 

l3 

q3 =1 
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Implementation (program # 5) 


1 The 64 QAM transmitted signal is given to both the I and Q channels through the For 
Loop 

2 91 (t) and 92 (t) are generated using QAM carr vi 

3 The mput signal is multiphed by the basis function and given to a Butterworth low 
pass filter of order 2 and sampling frequency set to lOOOOO (Adjustable based on 
carrier frequency) 

4 The output is mdexed at T/2 and stoied m indicators labeled RE and IM 

5 In sequence Structuie 2 this output of the LPF is given to multilevel case structure 
which will determme the position of the message symbol m the gray coded 
constellation 

6 The ii I 2 13 qi q 2 qs bits are determined usmg earlier defined decision rules 

7 The I and Q bits are mterleaved as ii qi 12 ,q 2 I 3 qa This was the 64 QAM symbol 
transmitted 

8 Output of the receiver is given to the "rec o/p " global variable and also to chart for 
plottmg 


Inphase channel 



quadrature channel 


Fig 2 9 64 QAM REC 
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23 4 NON COHERENT FSK 


The Quadrature realization of NCFSK has been implemented usmg correlators m 
the Inphase and Quadrature channels The correlator o/p are squared added and 
compared and based on this comparison a binary 0 or 1 output decision is taken[22] 

Implementation 

1 Generate The basis function V2/f cosrajt ■\j2IT sin co|t 
'v/2/T cosco2t and V2/T sm ©2^ 

2 Multiply the input signal with the above basis functions that is I and Q channels of ©i 
and ©2 

3 Give mdividual o/p to integrator Index integrator output at T 

4 Square the above obtained values and add polymorplncally 

5 Compare the two added values If the Inphase channel value is greater than the 
Quadrature channel value than the case structure gives a binary 0 as o/p otherwise a 
bmary 1 

NOTE The envelope detector realization of NCFSK is used in practical systems This 
detector has been implemented in Lab VIEW and is placed as an option in the select 
demodulator vi 


23 5 MINIMUM SHIFT KEYING 

The MSK signal is defined as 

s(t) = ai (t) cos — cos2Tcfct + ag (t) sin 

v2T; 


^ Tit ^ 


UTj 


sm 27ifct 


The MSK received signal is multiplied by the respective mphase and quadrature earners 
X(t) and Y(t) followed by the mtegrate and dump circuits The mtegration mterval is of 
2T seconds 
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s(l) = aj(t)cos 




X(t) = CO 




Ttt 




k2Tj 
COS27llct 


cos 27 tfct + aQ(t)sm 




2Tj 


sm 27ifct 


y 


Low frequency component of s(l) x X(t) equals 

Ttt 

ai (t) (1+cos— )/4 

The polarity of the sampler output determines the value of ai (t) bq (t) is deterimiied m a 
smular mannei[24] 



Fig 210MSKREC 


Implementation (program #6) 

1 Generate orthonormal basis function usmg MSK 1 Carr vi and MSK 0 Carr vi 

2 Add and subtract the two signals generated to get the basis functions X(t) and Y(t) 

3 Set counters 1 and k to mitial values 0 and 100 

4 Convert the incoming MSK signal from two dmiensional to a single dimension using 
Reshape Array function 

5 Take mput from local variable REC waveform over 2T i e 0 to 2T T to 3T, 2T to 4T 

and so on Correlate with X(t) and sample at 2T Give sampled values to INPHASE 
DECISION INDEX ^ 

6 After each mput of 2T mcrement L by 100 sample pomts which signify one T 
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7 For the Quradratiire channel take input over 2T from T to 3T 2T to 4T 3T to 5T and 
so on Multiply with Y(t) and pass to mtegrator Sample at 2T and give sampled 
values to QUADRATURE DECISION 

8 Increment k by 1 00 sample pomts which signify one T 

9 Give sampled values to case structure If INPHASE and QUARDRATURE indicator 
decisions are >0 then a binary 0 is detected If both are < 0 then also a binary 0 is 
detected For any othei combmation, a Bmary 1 is detected 

10 This piocessmg is done for the complete MSK signal produced by an array of 100 bits 
at the transmittei (considermg an mput sequence of 1 00 bits) 

1 1 The detected bits are given to chart for plottmg and to the global rec o/p for further 
processing 

NOTE The STOP button m all the modular VI s is to return to the Mam Menu 

24 LINE DECODERS 

The next stage after demodulation is to reconvert the bmary sequence to the origmal 
sequence by using appiopriate decoder operation A menu driven option for the selection 
of the decoder prompts the user to select the decoder An error message is displayed if 
the wrong decoder selection is made Decoders do the mverse operation of what the 
coders do The followmg decoders have been developed 

1 Biphase decodei 

2 Alternate Mark Inversion decoder 

3 3Bmary 4 Bmary decoder 

4 4 Bmary 3 Ternary decoder 

5 Bipolar Six Zero Substitution decoder 

The detailed description of each decoder is present m the VI information of mdividual 
programs However one decoder implementation m labVIEW is explamed below 
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241 AMI DECODER 



Fig2 11 AMI DECODER 

Ck =ak+ck-l 

bk = ck © Ck_i 

bk =(ak-f-ck_l)®ck_i 

With various combinations of a^ and cj^-i above, The following truth table is made 


ak 

bk 

4" 

1 

0 

0 


1 


Truth Table 2 7 (AMI decoder) 

IMPLIMENTATION (program # 1 ) 

1 Reinitialize all input variable to default values 

2 Using sequence structure 1 convert all mconimg bits to 0 or 1 using the truth table 
for AMI decoder 

3 We still have to extract the first bit Use wlule loop m sequence structure 2 (for 
the condition, Do While I = 0 and input bit is not equal to zero) 

4 If the input bit is = 1 then decode as 1 else decode as 0 
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5 


The previously decoded bits and this first bit are combined m the build array 
function to form the decoded AMI bit sequence 
6 Give to chart for plotting and to global variable "decoded o/p ' 

2 4 2 The next stage after decoding is the reconversion of the decoded bits to discrete 
voltage levels and to recover the input waveform, which was quantised modulated and 
transmitted The various options available for this operation are 

1 16 level linear decoder 

2 128 level linear decoder 

3 A Law decoder 

4 p. Law decoder 

The algorithm for A Law and p Law have been described below along with a detailed 
implementation of 128 level linear decoder 

A LAW ALGORITHM 

Discrete voltage level = (2Q+1) S=0 

= 2^(Q + 16 5) S=12 7 

where S is the segment and Q is the quantisation bracket Each 8 bit is mapped back to 
the corresponding voltage level S ranges from 0 to 7 and Q ranges from 0 to 1 5 (Refer 
to Table 2 8) 

u LAW ALGORITHM 

Discrete voltage level = (2Q+33) (2^) 33 

where S is the segment and Q is the Quantisation bracket Here too each 8 bit is mapped 
back to the correspondmg voltage level (Refer to Table 2 7) 
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128 Level Linear Recovery Implementation (Program # 8 ) 

1 The decoder o/p global output is converted to Boolean using For Loop and Case 
statement 

2 Convert one dimensional array of Boolean into a two dmiensional array (100*7) 
where seven bits sigmfy one voltage level 100 sample points have to be extracted 

3 Divide the voltage range by 128 and multiply with array (2° 2' 2^ 2^ 2 ^ 2^ 2^) 
which IS geneiated by usmg the For Loop and scale by power of 2 function Call this 
the multiplied array 

4 If first bit is 0 it signifies +ve voltage Use case statement for distinguishing between 
+ve and ve voltage +ve voltage being processed m case true and -ve voltage in case 
false 

5 Usmg shift registers multiply mput bits with multiplied array components 
polymorphically and add the one dunensional array components 

6 The resultant of the summation gives the voltage level pertaining to the 7 mputs bits 

7 Repeat process for all 100 one dimensional arrays of 7 bits and give the output voltage 
array to chart usmg the first bit for polarity information 

8 The mput signal which was mitially selected is replicated on chart with the 
quantisation error (The output signal is still digital It has been mterpolated and 
plotted on the graph ) 

Note For 16 level decoder same procedure is followed usmg 4 bits at one time 
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Figure 2 1 Function Generator 
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Binaiy No 

2 ® Compliment 

Sign Mag 

0000 

0 

+0 

0 

0 625 

0001 

1 

+1 

0 625 

125 

0010 

2 

+2 

1 25 

1 875 

0011 

3 

+3 

1 875 

2 5 

0100 

4 

+4 

2 5 

3 125 

0101 

5 

+5 

3 125 

3 75 

0110 

6 

-1-6 

3 75 

4 375 

0111 

7 (MAX) 

+7 

4 375 

5 V 

1000 

8(M1N) 


5 

4 375 

1001 

7 

1 

4 375 

3 75 

1010 

6 

2 

3 75 

3 125 

1011 

5 

3 

3 125 

25 

1100 

4 

4 

-2 5 

1 87 

1101 

3 

5 

-1 87 

1 25 

1110 

2 

6 

-1 25 

625 

nil 

1 

7 

625 

OV 


Table 2 1 16 level quantisation table 



Sc r c word 

Code word 

1 

000 

0101 

0101 

2 

001 

1001 

1001 

3 

010 

0001 

1110 

4 

on 

0010 

1101 

5 

100 

1000 

0111 

6 

101 

0100 

101 1 

7 

no 

0110 

0110 

8 

in 

1010 

1010 


Table 2 3 3B4B encoding table 
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Binary Inpul Block 

Ternary Output Block I 



Mode A 

Mode B 

Block Dipiial Sum 

0000 

+0 


0 

0001 

-+0 

-+0 

0 

0010 

0-+ 

0 + 

0 

0011 

+-0 

+.0 

0 

0100 

++0 

0 

±2 

0101 

0++ 

0 

±2 

0110 

+0+ 

0 

+2 

0111 

4-f+ 

i ** 

±3 

1000 

-H- 

+ 

±1 

1001 

-H- 

+ 

±1 

1010 

+-+ 

H- 

±1 

1011 

+00 

00 

±1 

1100 

OfO 

0 0 

±1 

1101 

0O+ 

00 

±1 

IhO 

0+- 

0+ 

0 

nil 

-Of 

0+ 

0 

T ible 2 4 4 Bimry 3 Ternary encoding table 
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Segment Code S 


Qua \ zal 
Code O 


Q anlizolion 
endpo nis 


000 

001 

010 

on 

100 

101 

no 

1 1 1 

0 

31 

95 

223 

479 

991 

2015 

4063 

1 

3 

103 

239 

511 

1055 

2143 

4319 

3 

39 

111 

255 

543 

1119 

2271 

4575 

5 

43 

119 

271 

575 

1183 

2399 

4831 

7 

47 

127 

287 

607 

1247 

2527 

5087 

9 

51 

135 

303 

639 

1311 

265d 

5343 

11 

55 

143 

319 

671 

1375 

.l733 

5 99 

13 

59 

151 

335 

703 

1439 

2911 

5855 

15 

63 

159 

351 

735 

1503 

3039 

61 1 1 

17 

67 

167 

367 

767 

1567 

3167 

6367 

19 

71 

175 

383 

799 

1631 

32 D 

6623 

21 

75 

183 

99 

831 

1695 

3423 

6879 

23 

79 

191 

415 

863 

1759 

3551 

7135 

2o 

83 

199 

431 

895 

1823 

3679 

7391 

27 

c7 

207 

447 

927 

1887 

3307 

7647 

29 

91 

215 

453 

959 

1 ol 

39 3 

7903 

31 

°5 

223 

479 

991 

2015 

4063 

8159 


0000 

0 

0001 

t 

0010 

2 

0011 

3 

0100 

4 

0101 

5 

0110 

6 

0111 

7 

1000 

c 

1001 

1O10 

10 

1011 

1 1 

1100 

12 

1101 

13 

1110 

14 

1111 

1 D 


Table 2 2 Law encoding table 


QiniMuHion 

ondpoinis 


Sr /nicnl Codr 


000 

001 

010 

Oil 

100 

101 

110 

111 

0 

32 

G4 

128 

2u6 

512 

1024 

2048 

2 

34 

68 

136 

272 

544 

1088 

2176 

4 

36 

72 

144 

208 

576 

1152 

2304 

G 

38 

7G 

1 j2 

304 

GOO 

1215 

2432 

S 

40 

80 

160 

320 

640 

1280 

2u60 

10 

42 

84 

168 

336 

672 

1344 

2688 

12 

44 

88 

176 

3j2 

704 

1408 

2816 

14 

46 

92 

184 

368 

736 

1472 

2944 

16 

48 

96 

192 

384 

768 

1536 

3072 

1 

uO 

100 

200 

400 

800 

1600 

3200 

•’C 

j2 

104 

203 

416 

832 

1664 

3328 

I’’ 

u4 

108 

216 

432 

064 

1728 

34 j6 

2 


112 

224 

448 

896 

1792 

3584 

2 

••o 

116 

232 

404 

920 

10 6 

3717 


( 

120 

240 

480 

JGO 

1920 

3040 


C2 

121 

240 

490 

992 

1994 

3 )08 

32 

04 

128 

2j6 

u12 

1024 

2018 

4096 


Qu nl z iU( 1 
Codd 


0000 0 
0001 1 
0010 2 
0011 3 

0100 4 

0101 5 

0110 6 
0111 7 

1000 
1001 9 

1010 10 
1011 11 
1100 12 
1 101 1 i 

1110 11 
1111 1u 


Table 2 8 A Law encoding table 
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CHAPTER 3 


BASEBAND LINEAR EQTTAl J/RRS 

111 the previous chapter we saw the LabVIEW implementation of various passband 
modulation & demodulation schemes In this chapter we simulate ISI usmg different 
channel characteristics and pass this ISI effected signal through a Zero Forcmg equalizer 
then a preset type equalizer which uses the LMS algorithm to find the optimal tap gam 
coefficients and lastly through a hybrid equalizer one that is a combmation of a preset 
and an adaptive equalizer The mcrease m the noise margms and the effect on detection of 
the mput data through these equalizers have been highhghted through these simulation 
processes The detailed implementation of the preset equalizer is explamed with its 
LabVIEW block diagram 

One of the main factors which contribute to Inlersymbol mterference is the 
channel characteiistic In Pulse Amplitude Modulation the Intersymbol mterference (ISI) 
results firom Imear amphtude distortion m the channels that broadens the pulse and causes 
them to mterfere with one another The three channel impulse response characteristics 
used m sunulation are given m Figure 3 1 The first phase of simulation deals with the 
effect of the channel impulse response and the resultant ISI so produced on the 
probability of error of detection The eye diagram and the BER plots which constitute 
the second phase are discussed m chapter 4 

3 1 Zero Forcing Linear Equalizers 

The total number of taps in the tapped delay line equalizer implemented is 
(2N+ 1 ) The impulse response of the Tapped delay Ime is 
N 

X WK5(f - kt) S(t) = Dirac delta function 

h(t) = 

and = weight at k"' tap 
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If this tapped delay filter is cascaded with a system with impulse response c(t) 
then the overall response is [13] 

p(t) = c(t) * h(t) where * stands for convolution 
= c(t)* 'Zw,S(t-kT) 

k A 

k- N 

at t = nT samplmg instants 

P(nT)= Y^Ai^-k)T) 

k—N 

For minimizmg ISI 

p/ -rN _ 1 1 11=0 

^ [0 n^^O n=±12 +N 



Figure 3 1 Channel Impulse Response 
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If Cn = c(nT) then the function of the Tapped delay line can be described as 



- c 

— C — 1 



"o' 

- 

- - 

- - 

^ +. 


0 


- Co 

- c 

w 

= 

1 


- - 

- - 

w, 


0 

p2/ 

- C„ 

- 



0 


The larger the numbers of taps the lesser is the distortion 

Implementation 

1 Geneiate landonily a binary sequence of suitable length and convert to pulses with 
time duration T 

2 Select any discrete channel from the three characteristics m Figure 3 1 and 
convolve each bmary pulse with the channel Call this h(t) 

3 Generate response of the matched filter which matches the channel response Call 
this h( t) 

4 Convolve h(t) and h( t) Repeat this for the complete sequence of input bits 

5 Consideimg only 9 taps on the equalizer sample each h(t)*h( t) output at its peak 
and at 9T 3T 2T T T 2T 3T 9T Form a 2 dimensional array with 
these discrete samples for all mput bits The basic assumption is that the amplitude 
at other sampling mstants is negligible 

6 Shift each one dimensional array by one T as we progress down the two 
dimensional array That is 2"' Row shift by 1 T 3^' row shift by 2 T and so on 

7 Add the discrete values m each column and this is the combmed Imear effect of 
ISI What has been described in the above steps is that each pulse is treated 
mdividually and these individually effected pulses are Imearly added to simulate the 

overall ISI effect 

8 The samples C(nT) are then used to solve the simultaneous equations as described 
m (I) to determine w, w 3 w. wjw, Smee N=4 so the dimension of the C 

matrix is 9 x 9 
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9 


On knowing the coefficients of the Tapped delay filter the output samples are 
determined by 

N 

Po[(k+N)T]= i:CnPr[(k-N)T] 
n=-N 

where Po (t) is the output sample and Pr(t) is the input to the tapped delay filter 

The assumption here is that the pulse broadened due to the channel will have zero 

mterfeience at sampling mstants 4T before and 4T after the peak 

10 Each stage of the above implementation is displayed on graphs generated m the 
labVIEW program 

PLOT 1 is the mput bmary symbols 
PLOT 2 IS the Channel and transmitted pulse convolution 
PLOT 3 IS the Matched filter impulse response 
PLOT 4 IS the Convolution of h(t) and h( t) 

PLOT 5 IS the detected bits 

1 1 Fmally the number of errors between the detected bits and the transmitted bits are 
determined 

12 The detection errors with all three channels with varymg number of mput bits are 
discussed m chapter 4 Provision for adding AWGN samples to the transmitted 
data symbol signal is also available 

3 2 PRESET EQUALISERS 

In a preset equalizer the error components are measured by transmittmg test pulses 
(PN sequence of maximal length 5 used m the implementation) through the system and 
mspectmg the output of the transversal filter at the samplmg mstants The equalizer is 
adjusted prior to or durmg breaks m data transmission by these predetermmed test pulses 
The Zero Forcmg equalization assumes that the imtial distortion is less than 100 percent 
although higher values could be encountered m practice The second approach the rms 
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equalization where the tap gams must be adjusted such that the cross correlation between 
the error voltage e,, and the mput samples x„ at delays of N N is zero 


The mean square distortion of the output response is defined as 


2 

s = — -T h 2 
2 n 


ho" n 

Constrammg ho usmg Lagrange multipher we get [21] 

00 

Y^2h^x^_j=Xj whe}eJ = -N N 


since 


J -N CO 

N CO 

,=Xj whcfcb^ = Jfc 


A- N 


(I) 


(II) 


The best tap gams m the mean square sense are therefore obtamed by solvmg 
simultaneously (2N+1) linear equations (II) 

Let e,, denote the error at sample n present m the output pulse h(t) as compared to 
the ideal unity height pulse 

6n = hn n 0 (III) 

€o = ho 1 elsewhere 

Then from (I) and (III ) we get 

00 

;^€nXn_T =0 where j=-N N 
n=-oo 

This proves that the cross correlation between error signal and input pulse should be 
zero at each sample point withm the range of the equalizer 

Implementation 

1 Generate PN sequence with m = 5 and mitial state of shift registers as 1 0 0 0 0 
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2 Generate baseband pulses of duration T with Amphtude +1 for a binary 1 and 1 for 0 

3 Select channel impulse response A (see figure 3 1) Convolve input pulse with channel 
impulse 1 espouse call tins h(t) 

4 Generate matched filter i espouse with the channel impulse response Convolve with 
h(t) 

5 Preset equalizer being a 7 tap equalizer (In simulation) the discrete amplitude at 3T 
2T T T 2T 3T IS stored m array form for each mcommg binary pulse Then each 
row of this 2D array is shifted and addition in each column is done The resultant are 
the ISI effected samples at mstants nT (Each pulse effects 3 precedmg and 3 
succeeding pulses ) 

6 The output of the matched filter is given to the tapped delay filter wheie the weights 
are calculated by the LMS ( least mean square) algorithm[18] The weights determined 
by the LMS algorithm makes the cross correlation between the error signal and the 
input pulse go to zero 



Figure 3 2 LMS Algorithm 


Wk+1 = + 2pskXk 

where p is the gam constant that regulates the speed and stability of adaptation Wk is 
the tap gam vector and Xj^ is the mput vector to the equalizer filter The LMS algorithm 
runs iteratively for a selected value of p till such time the error oscillates around zero The 
weight vector converges to the optimum vector solution 
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7 Generate a suitable length binary sequence call this the data sequence FoUowmg step 
2 through 5 generate the resultant disci ete ISI effected samples which will be the input 
to the Tiansversal filter 

8 With the optimum weights determined by the trainmg sequence carry out operation 

yn=xjw„ 

wheie = [xn+3 ^n+l Xn-i 

fV =[w W W3] 

yn=Lyn+3 yn yn- 3 ] 

9 The decision is based on the yn The threshold is set to zero The detected bits are 
compared to the input data bits foi errors Gaussian noise can also be added to the 
data bits 

1 0 Output graphs at each stage as mentioned earher can be viewed at each stage in the 
labVIEW program 

3 3 ADAPTIVE EQUALIZATION 

In adaptive equalization the error voltage €„ are continuously estimated during the 
course of normal data transmission and corrections to the equalizer coefficients are 
effected as required [21] Such closed loop system has an inherent accuracy advantage 
over the preset type of equalizer In random data transmission the transmitted sequence is 
not available but the estimate at the output of the threshold detector is used as if it were 
the correct transmitted sequence {an} The error voltage is given by 

Sn ~ yn 

where yn is the equalizer output at time nT and an is the estimate after the threshold 
detector 


50 




N 


X 


on AY T 



DELAY T 


DELAY T 



DELAY T 







1 





T 1 

T 

r- 

5 Wo ^ 

^ 1 

5 

5 0 

+ 1 

yn 


— - — - - t 

^ e 


Figure 3 3 Elements of an Adaptive Filter 

Implementation of Decision Directed Equalizer 

1 Step 1 Ihiough 6 in the pieset equalizer implementation are same for the decision 
directed mode equalizer Once the optimal tap coefScients are determined generate a 
suitable length of randomly generated bmary numbers usmg the iterative For Loop 
and the random no generating function 

2 Each pulse is transmitted over the channel and the resultant mput to the equalizer is 

C(t)*h(t)*h( t) Where * stands for convolution 

Y 1 = X,/ Wn where Xn are the mput to the equalizer Wn is the tap gam coefficient 

vector and Yn is the response of the equahzer 


3 The threshold detector (case statement) gives output a„ 1 or 1 The difference an yn 
is determined and the tap weights are readjusted usmg this value of en The modified 
weights would result m a lower error en 

4 With each detection of the transmitted pulse the correction to the weights is done 

5 The impulse response of the channel was varied at a slow pace by a magmtude of 10 ^ 
at each discrete point of it s response with each transmitted bit In one case the 
channel was tracked and m the other left untracked With varymg noise variance it 
was seen that the data transmitted through the non tracked channel gave 5 times more 
errors than the tracked channel where the optimal weights were changmg dynamically 


lENTRAL UBRAr 
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This program mcludes plots of this varying channel and the error indicators of the 
tracked as well as the untracked channel The basic stage wise plots as before are 
present as well 

The Pieset channel vi (program # 9) 

1 Sequence structure 0 is used for initializing all mput variables to default values 

2 Sequence 1 generates the PN sequence and converts the generated symbols to polar 1 
and -1 using baseband vi AWGN Noise samples are generated which are added to 
the data samples 

3 Sequence 2 convolves the input pulse and the channel response and further convolves 
the matched filter response with the resultant of the above convolution 

4 Sequence 3 shills the output at the matched filter by the corresponding T and then 
tiansposts this 2 dimensional aitay and adds m a row wise fashion It generates the 

input to the equalizer liltei ^ 

5 Sequence 4 first runs the LMS algorithm and determmes the optimal weights The 
mput vector to the equalizer filter and these weights are added and a detection 
decision is taken whethei a bmaiy 1 was transmitted or a binaiy 0 

6 Sequence 5 gener ites the data bits using the random number generation function The , p- 

length of the binary sequence is specified by the no of bits control Step 2 through 

4 repeat for the data sequence and the output of the matched filter is given to the 

equalizer 

7 Sequence 6 multiplies this input vector with the optimal weights determmed earher 

8 The output of the sum array element function is given to a case statement where the 

binary sequence is reconstructed 

9 In case of the adaptive equalizer vi the difference between a„ and y„ is given as the 
error value to correct the optimal weights The next set of mput to the 2N+1 tap 
equalizer filter use these changed gam coefficients This is a closed loop process where 

after each detection of a bit the weights are adjusted 
The various processing stages of the preset vi are depicted in the flow chart placed at the 

appendix 


52 


CHAPTER 4 


PARAMETER DETERMINATION AND SIMULATION 

RESULTS 


In the pievious two chapters we have discussed the simulation of various 
passband moduhtion schemes and equalizers m the baseband case In this chapter 
Section 4 1 we discuss the various parameters to be determmed such as the 
amplitude spectium, power spectrum of the mput signal and the modulated carrier 
signal the lesponse of the filters used m the detectors the Gaussian pdf of the noise 
samples and the quantisation parameters The plots of these above parameters 
obtained through LabVIEW implementation are also discussed The BCR plots of 
some moduhtion schemes are discussed in Section 4 2 In Section 4 3 the output 
signal gcncia'ed by QPSK OQPSK and MSK schemes are compared and a 
compiiison of the eiiors geneiated by bandlimited OQPSK and QPSK modulation 
with vaiying noise van incc is discussed The eye diagram of the equalizer output and 
the effect of channel characteristics are discussed m Section 4 4 Section 4 5 deals 
with the regeneration of an externally generated signal usmg BPSK modulation 
scheme, on a CRO 

41 PARAMETERS 

In the passband modulation demodulation capsule following parameters of the 

output ptocesses aie determined 

I ) Amplitude and power spectrum of the input signal 

II ) Amplitude and power spectrum of the modulated carrier signal 

III) Quantisation parameters 

IV ) Filter response 

V ) PDF of AWGN noise samples 

I 

4 1 1 Amplitude and Power S pectrum 

The amplitude and power spectrum is computed as 

FFT(Singal) 

(N) 
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where N is the number of samples points in the signal The program then converts the 
amplitude to single sided rms magnitude spectra Plot I shows the amplitude spectrum 


of a 10 kHz signal 



1 he power spectrum is computed as 

rrT * (Singal)xFFT(signal) 

(Nh 

where N is the number of points in the signal and * denotes the complex conjugate 
The piogram then converts the power spectrum into a single sided power spectrum 
Plot II (a) shows the power spectrum of a MSK transmitted signal at 80000 Hz 

dBVrmt | 02 ' 

ao 

40 

60 
75 

80000 0 90000 0 100000 0 110000 0 120000 0 130000 0 140000 0 1 50000 (1 



. , MSK « O D U LAT E 0 C A B R IE R I 

(a) 



(c) 


PLOT II 









The plot II (b) and (c) above show the power spectrum of a QPSK signal at lOOOOO 
Hz and 16 QAM modulated carrier at 10000 Hz 

Note Both the power spectrum plots use 1500 samples 

In all the plots above the units and the scales (logarithmic / linear) can be selected as 
per the options available on the front panel 

412 Quantisation parameters 

The following pnnmeters of the linear quantizers used are determined m the 
quantiser parameter vi and are displayed on its front panel using slide indicators and 
formula nodes 

I ) a (spacing) 

II ) Mean square error (a^/1 2) 

III ) Root mean square error (a/ Vn) 

where M is the level of quantisation 
V ) Signal to noise ratio (peak/rms) V/(a/-\/l2 ) 

4 1 3 Filter response 

This program graphically shows the transfer function of a butterworth filter 
used as a band pass or a low pass filter The order sampling rate and cut off 
frequencies (low and high) have to be specified as inputs parameters to the programs 
These parameters are set to some default values The functioning of the program is as 
follows 

1 Generate 1024 samples of impulse pattern 

2 Pass through desired fi Iter (Low pass or Band pass) 

3 FTT of filtered output is converted to polar form 

4 The magnitude gives the impulse response plot and the angle gives the 
phase response plot 

Plot III shows the magnitude and phase spectrum of a bandpass butterworth filter The 
order and the cut off frequencies are set to some default values 


iv ) Mean signal power (aQ 


M ^-1 

12 
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4 1 4 Distribution and Density Function of AWGN samples 

This program determines the probability distribution of the AWGN samples 
using the histogram and integration functions available m LabVIEW The density 
function IS determined by taking the derivative of the distribution fiinction The 
minimum samples required to generate a bell shaped density function is 1000 Plot IV 
gives the pdf and the distribution of 1000 samples of AWGN The input samples are 
also plotted along with the pdf 
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4 2 PROBABILITY OF ERROR 

For digital communication system the relevant measure of performance is 

always related to the systems error producing behavior BER is the fractional 

number of errors in a transmitted sequence The conceptual framework within which 

this problem can be visualized is as follows The mput is defined as an ensemble of 

identical systems the input to each of which is a distmct joint realization of the signal 

and noise source The input to the decision device is this set of waveforms and the 

output IS a digital sequence which in general will contain some errors For the kth 

member of the ensemble the error rate is defined as 

Pw= hm nk(N)/N 
N->oo 

Not all Pk will be identical since these are realizations of the digital sequence which 
are not typical The setting within which P is defined is identical to the classical 
Bernoulli trails with probability P of success and (1 P) of failure 

"Monte Carlo is merely a label for the implementation of a sequence of 
Bernoulli Trials in which we count the number of errors and divide by the number of 
trails This method requires no assumptions about the input processes or the system 
BER for the assumed symbol can be written as 

P = (ueDo fv(u)do 

where fy is the pdf of the sampled symbol at sampling epoch and Do is the region of v 
corresponding to an error 

Define error indicator function as 



P=r jr(v)A(u)cio 

*>-00 

P=E[H(v)] 

The estimator P of the expectation is the sample mean we have 

lei 

where No is the number of elements and V, is the sequence of symbol space samples 
of the decision voltage lo is the integer set that contains all i 
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I 


If all symbol have the same probability of occurrence then 
F = n(N)IN 

where N are the total symbols & n is total errors 
As N->oo P -> P a s 


To generate the BER graphs for various modulation schemes the following guidelines 
and assumptions are used 

1 20000 data bits are generated from a random source for processing 

2 AWGN samples with varying seed are added to the input waveform generated 
from the data bits 

3 Averaging over 5 times at each input variance value is done to smoothen the 
BER curve 

{Due to limitations of processor larger set of data bits could not be 
considered} 

4 Eor the correlator receivers No/2 = cr^ where No/2 is the single sided noise 
density in watt/I Iz [23 ] 


4 2 1 Probability of Error BPSK 


PcBi SK = ^erfe 


Eb 

N, 




0 J 


Prictical Implementation 

BPSK Error vi (program # 10) runs the PSK Verific vi, which calculates the 
number of bits in error at varying input noise variance 


1 Randomly generate input signal, which is quantised by a 4 level linear quantiser 
producing 400 bits per input signal 

2 The input bit stream is given to a BPSK modulator which generates the carriers 
for the binary 0 and 1 

3 AWGN samples at desired (varying) variance is added to the carrier and this 
carrier plus noise is given to the receiver 

4 The correlator realization of the BPSK receiver detects these bits according to the 

polarity at the output of the threshold detector 
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5 The input bit sequence is compared to the detected bit and the number of bits in 
error are determined for each run The number of bits selected is between lOOOO 
to 20000 At each variance averaging is done over 5 runs 

6 The average signal power to average noise power (SNR) is dB is the x axis of the 
plot with the y axis being the bits in error to total bits The plot is in logarithmic 
scale 

7 The BER plots had to be extrapolated for higher values of SNR due to limitations 
of simulating a run with larger bit sequence The processing time for a sequence 
1 0000 bit long would take upto 3 to 4 hours 

8 The theoretical plot of Pe and Eb/No is plotted (plot#l) along with the practical 
plot using the expression 



BER PLOT 



PLOT 0 1 BP5K Et>''N 

N OF BITS TRANSMITTED 12000 
PROCESSINQ TIME 3H 


NOTE In plot #1 to plot # 4 plot 0 is the theoretical plot and plot 1 is the simulated 
plot 
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422 Probability of Error OPSK 


P 


QPSk ~ 2 


No 


QPSK Innsmits it twice the bit rate of a coherent binary PSK system for the same 
channel bandwidth 

The plots of BER Vs SNR obtained practically and theoretically are given in Plot#2 
The QPSK error vi runs the same way as BPSK error vi with QPSK modulator and 
receiver 


4 2 3 Probability of Error DPSK 


Pe 


1 f-E^ 

—exp 

2 UNoj 


E = 2Ei, 


No = power spectral density 


Pe 


1 

= -exp 


' -Eb 
, No 


\ 

y 


Eb = signal energy per bit 


The Simulated and the theoretical plot are given m Plot#3 The practical 
implementation has been done using DPSK Error vi using DPSK transmitter and 
receiver with AWGN samples 


BER PLOT BER plot 



,PI to 

r 

iPlotl 


dB 


PLOT *1 0P«K 

M or BITS TRAKIMITTED HOOD 
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PLOT P3 DPSK 
N OFBITSTRANSMITTED 


12000 



^ Probability of Error Non Coherent Frequency Shift Keying 

Plot#4 gives the piobability of error Vs SNR plot of noncoherent FSK found 
piacticilly and theoretically 

'Ihc theoretical plot is found using 


Pe =-^exp(-Eb/2No) 



4 2 5 Probability of Error 16-OAM 

The positions of the bits in the QAM symbols associated with each pomt m 
QAM constellation have an effect on the probability of them being m error Of the 
four bits in 16 QAM, for i] and qi (MSB), the distance from the demodulation 
decision boundary of each received phasor in the absence of noise is 3d for 50% of 
the time and d for the rest 12 and qz (LSB) are always at a distance d 16 QAM may 
be considered as having class Cl and class C2 subchannel where Cl has lower 

probability of error than C2 [19] 
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I or C2 subch uiiiel 


r d 1 °° 

Eo ~ 1 Od [ Square constellation with gray coding] 
tor Cl subch innel 



Combining P 20 and Piq 

Pav = (PlG + P2a)/2 

The theoretical plot has been based on the above expressions 

Passing bits tlirough the 16 QAM modulator AWGN channel and 
demodulator ch un and checking for errots at varymg noise variance we obtain the 
practical icsults Ihe seed of the AWGN sample generator has been varied in the 
generation of the gnussian samples Plot#5 shows the BER for 16 QAM 


NOTE The theoretical plot is the symbol error Vs Eo/No plot where Eo is the 
average energy and No is the spectral density of noise The low pass filter realization 
of the detector which is used is not the optimum detector realization 
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>PSK, QQPSKAND MSK SIGNAL GENERATION 


4 3 1 In this section we compare the QPSK signal OQPSK signal and the MSK 
signal, gcnciated at the lespectivc tiansmitters 

1 ho MSK sign il has no tiansitions of phase and it is continuous m phase 

I ho QPSK signal has a maximum phase transitions of 180° if both Inphase 
and Quardiatuic oomponont change phase simultaneously 

The OQPSK can have a maximum phase transition upto ±90° due to the offset 
by Tb m the inphase and quardrature channel 
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4 3 2 QPSK, OQPSK Through Bandpass Filter 

The histogram plot VI below give the comparison of the errors generated in 
the QPSK modulation scheme and the OQPSK modulation schemes by varying the 
gaussian noise viri nice and passing them through a bandpass filter Comparison with 
a run of 7000 bits and 5000 bits are shown m the plot It is observed that OQPSK 
modulation scheme generates fewer errors than the QPSK modulation scheme 

A 70o”o bif long input binary sequence is randomly generated (fl and f2 for the 
carriers is 100000 Hz High cut off frequency of the bandpass filter is 195000 Hz 
Low pass frequency cut off is 5000 Hz The order of the butterworth filter is 2 ) 
This IS given to the QPSK transmitter, which generates the QPSK signal This 
signal IS bandpass filtered and given to the detector AWGN noise samples are 
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1 ) 

2 ) 



added to the transmitted signal At the detector the coherent earner so e ener„i,H 
IS also passed through a bandpass filter which has the same high frequency cut off 
and low frequency cut off as the QPSK signal 

3) The detection errois at varying noise vaiiance are recorded m a spreadsheet format 
file 

4) The OQPSK signal generated is passed through a bandpass filter and detected at 
the detectoi Here too the coherent carrier so generated is passed through a band 
pass filter with the same parameters The detection errors at varymg noise 
variance are recorded m a spreadsheet format file 

5) The errors at various noise variances are compared for QPSK and OQPSK a 
histogiam plot is generated foi a direct comparison 

6) It lb observed that the OQPSK generates less errors then QPSK 
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4 4 EQUALIZERS AND EYE DIAGRAM 

4 41 Channel Impulse Response and Errors Due to ISI 

The channel impulse response plays an important role m the generation of ISI 
In some channels the ISI can be so severe that even in the absence of noise detection 
errors persist Figure 3 1 gives the three channel impulse response used m the 
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simulation The ISI pulse so constructed (Refer Implementation of Zero forcmg 
bquali/er) is given to the ZF equalizer usmg all tliree chaimel characteristics The 
eiiois detected with vaiying input bit size is recorded m ZFE spreadsheet file 
Table 2 8 bungs out the detection enors for various runs of the simulation for all the 
thiee chinncl impulse response 
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Table 2 9 Bit errors at various runs 


4 4 2 Eve Diagram Generation 

Refer to implement ition of adaptive equalizer in Chapter 3 The signal before 
the transversal liltcr and the signal after the transveisal filter are used to generate the 
eye patterns be lore and iftei equalization The mterpolated signal has been 
fragmented into signal over T sec These T sec signals are supeiimposed (stacked) 
one over the othei I he input signal ovei 100 T (100 bits generated) when stacked 
generates the eye diagram before equalization The eye diagram after the equalizer is 
generated usmg the output signal from the equalizer (Interpolated) and again 
superimposing (stacking) 1 sections of the signal It is obseived from the eye 
diagiams simulated tint theic is an impiovement in the noise margin in the signal 
after the transversal filter Plot VII shows the eye diagram before and after 
equalization The data for eonstructmg the eye patterns has been extracted from the 
Preset equalizer vi 
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EYE PATTERN BEFORE EQUALIZATION 






PLOTVn 

4 4 4 Probability of Error For PAM With ISI 

II a limo hmituci signal is limited to N L ary symbols the ISI can only assume 
r ^ possible values However il the impulse response x(t) is time limited to k symbols 
belore md viler to then sequences need to be considered to evaluate the 
probability ol crroi 1 he probability of error can be determined by averagmg over all 

sequences (21] 

1 2k 

pYiPO) 

P„ IS the probabtllly of error for the i* sequence and P(i) is the probability of 

occiirrenec for the i'*' sequence 
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where D(i) is the distortion over ith sequence 
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Pl = P [\D(i) + ; 7 o| > V d] where rjo is the noise sample and Xod is 

the distance from threshold 
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where the noise is gaussian with zero mean and variance cr^ [23 ] 

D(i) has been evaluated for each input sequence to the adaptive equalizer in the 
adaptive equalize! vi The threshold distance Xod has been taken as 1 The 
theoretical Pe has been calculated by averaging over all input sequences The signal 
amplitude is 1 the detection errois at varying noise variance have been evaluated for 
a 1 000 bit binary input sequence 

Simulated graphs of the practical implementation and theoretical implementation are 
given m plot VIII 
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4 5 REGENERATION OF EXTERNAL SIGNAL AFTER 
PROCESSING 


A sinusoid'll sign'll nt IkHz fiequency is taken fiom a flinction generator 
through DAQ interfice card and DAQ function generator The DAQ function in 
LabVIEW C Scan is used to scan 15 cycles 100 sample points each at a sampling 
rate of 1 0000 from channel 1 of device 1 (function generator) This input signal is 
quantised and BPSK modulated The BPSK receiver and decoder regenerate this 
IKHz signal This reconstructed signal through the DAQ interface using analogue 
out continuous generition flinction is given to the DSS 5040 oscilloscope The 
output signal on CRO is a IKIIz signal, the quantisation steps at discrete levels are 
visible on the CRO scope in the digital storage mode 
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CHAPTF.R 5 


CONCLUSION 

rht sirnul'ition system which has been described allows interactive digital 
simiilition ol ‘iibitriiy point to point digital communication system The software has 
been stiuctuicd to allow flexible interactive use by even the most casual usei The 
iTi ijol ipplic ition of this communication simulator package will be m support of new 
development efforts m modulation/coding concepts which are emerging anew m the 
communication woild The various system parameters the graphical representation 
of the wivefoims at any stage in the link which can be displayed with ease give the 
usei a bettei undeist mdmg of i practical system Spectrum analyzers or CRO s are 
not needed is one can view these wavelorms oi their spectrum on the LabVIEW 
displ lys 1 he gener il philosophy in developing the software has been that it should 
be used is m idjunet ind not i substitute of the system 

1 abVll W with Its debugging features the graphical support and its mterface 
for dati acquisition mikes it a vciy suitable choice foi using it for simulation The 
other simul ition paekiges like simuhnk and higher level languages like FORTRAN 
could also h ive been used but the basic advantage of a graphical mterface and data 
acquisition would not be picsent The various data operations plots/charts along with 
the vmous lunetions and library ol Vis in LabVIEW makes the simulation job much 
more ittraetive and easiet Simulation in LabVIEW takes one closer to the actual 
hardwaie implement Uion ol the system in hardware 

Ihe simulation results discussed in chapter 4 do conciu with the results 
obtained m theory IK) H 14 22] The BER plots would have been closer 
approximations had a largei data sequence been considered and averaging over a 
larger number of runs been carried out The processor limitation can be abstamed by 
doing these simulation runs over a mam frame 

Another major advantage of the simulation package is the ease with which it 
can be updated to incorporate enhanced capabilities or to add new modulation 
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schemes Modules of these new modulation schemes can be developed and easily be 
accommodated The ability to conveniently expand system capabilities will ensure 
that this system will provide a usefiil tool in support of these efforts 


SUGGESTIONS FOR FURTHER WORK 

In Section 4 5 it has been indicated tliat leal time signal processmg is possible 
using the DAQ interface If a dedicated host computmg machine with a faster 
processor were used it would be possible to sunulate a real tune system A 
transmittei and reccivei model could be developed using some modulation format to 
simulitc and itcovci the tiansmitted data A real time baseband communication 
system could be developed wheie the tiansmitter could be on one PC and the leceiver 
on inothei 

I he modulation toimat which have been developed are the basic ones modem 
niodulition techniques used like TCM OFDM could be incorporated Multiple 
access techniques like fDMA FDMA can also be included as simulation hinctions 
Carrier and bit synchrony Uion can be included in the simulation feature by designing 
a 2nd ordci PLL module Random error correcting codes like convolution codes 
BCn codes and othei such schemes can also be provided 

The channel assumed is a Gaussian channel A practical channel would consist 
ol a lading multipath (dispcisivc) channel along with the AWGN channel A fadmg 
channel could be mcorpoi itcd in the simulation The distortions due to narrow band 
filtering c in ilso be included as a function of the transmitter In the baseband 
schemes raised cosine filters at the transmitter and receiver could be implemented 
The equalizers which have been designed for 9 taps could be extended to 65 taps 
(gam coelficicnt m pi ictical systems) and more discrete channel model may be 

mcorpoi ated m simulation 
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